SHARP BOUNDS FOR HARMONIC NUMBERS 
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Abstract. In the paper, we first survey some results on inequalities for bound- 
ing harmonic numbers or Euler-Mascheroni constant, and then we establish a 
new sharp double inequality for bounding harmonic numbers as follows: For 
n £ N, the double inequality 

<H(n) — Inn 7 < 



12n 2 + 2(7 - 127)/(27 - 1) ~ 2n 12n 2 + 6/5 

is valid, with equality in the left-hand side only when n = 1, where the scalars 
2 ^27~i~^ an< ^ § are best possible. 



The series 



1. Introduction 



1 + I + I + - + I + (1) 



is called harmonic series. The n-th harmonic number H{n) for n G N, the sum of 
the first n terms of the harmonic series, may be given analytically by 



I 

H{n) = Y j - = 7 + V(n + l), (2) 

»=i 

see [1, p. 258, 6.3.2], where 7 = 0.57721566- • • is Euler-Mascheroni constant and 
ip(x) denotes the psi function, the logarithmic derivative j^y of the classical Euler 
gamma function T(x) which may be defined by 

Y(x) = / e^e^dt, x>0. (3) 







In [17], the so-called Franel's inequality in literature was given by 

<if(n)-lnn-7<^-, n G N. (4) 
2n 8rr 2n 



In [11, pp. 105-106], by considering 

.1 - 



In = 



1/r. 



X 



dx = lnn — H(n) (5) 
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and < /„ < |, where [t] denotes the largest integer less tan or equal to t, it was 
established that 

^ <H{n)-Ynn< 1, n £ N. (6) 
In [11, pp. 128-129, Problem 65], it was verified that 
1 ™ 1 1 

-ln(2n+l)<^— - T <l + -ln(2n-l), n £ N. (7) 
fe=i 

In [27], it was obtained that 

L — < ff(n)-lnn-7 < n e N. (8) 

1 2n 



2(n+l) 
In [8] , it was proved that 



! <tf(n)-ln(n+M - 7 < , n £ N. (9) 



24(n + 1) 2 v 7 V 2/ 24n 2 ' 
In [21], the following problems were proposed: 

(1) Prove that for every positive integer n we have 

^<^)-lnn- 7 <^_. (10) 

(2) Show that | can be replaced by a slightly smaller number, but | cannot 
be replaced by a slightly larger number. 

In [10], by using 

^ = ln ^ + i-ii + i2^ <") 

for < s n < 1 , these problems were answered affirmatively. The editorial comment 
in [10] said that the number | in (10) can be replaced by 2 ^~* and equality holds 
only when n = 1. This means that 

1 1 

<if(n)-lnn-7< - — — , n £ N. (12) 



2n + ^ - 2 " 2n + 3 

This double inequality was recovered and sharpened in [6, 7] and [18, Theorem 2]. 

In [26] , basing on an improved Euler-Maclaurin summation formula, some general 
inequalities for the n-th harmonic number H{n) are established, including recovery 
of the inequality (10). 

In [25], the problems above-mentioned was solved once again by employing 

TT< \ 1 J- a. 1 1 B * f°° ^M. A fl* 

^ (n)=lnn + 7+ 2^-2g^-/ n (13) 



J n 



and 

X 2 q aX< 2qn 2 q ' 

where n and q are positive integers, Bi{x) are Bernoulli polynomials and B 2i = 
B 2 i(0) denote Bernoulli numbers for i £ N. For definitions of Bi{x) and B 2il please 
refer to [1, p. 804]. 

In [23], the inequality (10) was verified again by calculus. 
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In [13], by utilizing Eulcr-Maclaurin summation formula, the following general 
result was obtained: 

Hin) - L, n + 7 + - - — 2 + - g ( , _ + o(^) . (15) 

See also [15, p. 77]. In fact, this is equivalent to the formula in [1, p. 259, 6.3.18]. 
In [22], by considering the decreasing monotonicity of the sequence 

Xn = — rivOTi ~ 2n ' (16) 

it was shown that the best constants a and b such that 

oo 1 



1 

< 



2n + a 



< 2^6 < 17) 



for n > 1 arc a = 1 J\ n2 — 2 and 6=1. 

In [4, Theorem 2.8] and [19], alternative sharp bounds for H(n) were presented: 
For n e N, 

1 + ]n(y/e - l) - In(e 1 /(" +1 ) - l) < H(n) < 7 - ln(e 1 /( n+1 ) - l). (18) 

The constants 1 + ln(y / e — l) and 7 in (18) are the best possible. This improves 
the result in [3, pp. 386-387]. 
In [20], it was established that 



In + + 7 < H(n) < ln(n + e 1 7 - l) + 7, n £ N. 



(19) 



In [5], it was obtained that 



1 <H{n)-ln(n+l ) - 7 < — % ..„ (20) 



24{n + 1/2 V /6[1 - 7 - ln(3/2)] } " V 2 / 24(n + 1/2) 

for n £ N, where the constants 



2^611-7-^(3/2)] 

and ^ are the best possible. 

For more information on estimates of harmonic numbers H(n), please refer to 
[9, 24], [14, pp. 68-86], [15, pp. 75-79] and closely-related references therein. 

The aim of this paper is to establish a double inequality for bounding harmonic 
numbers, which is sharp and refines those inequalities above-mentioned. 

Theorem 1. For n £ N, the double inequality 

<g(n)-lnn-^--7<- — / — (21) 



12n 2 + 2(7- 12 7 )/(2 7 - 1) ~ w 2n ' 12n 2 + 6/5 

is valid, with equality in the left-hand side of (21) only when n = 1, where the 
scalars 2 ^~^ 7 - 1 and | in (21) are the best possible. 

Remark 1. When n > 2, the double inequality (21) refines (20) and those mentioned 
before it. 
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2. Proof of Theorem 1 

We now are in a position to prove Theorem 1. 
Let 

= In. + l/J-^ + l) - (22) 
for x € (0, oo). An easy computation gives 

4xV(a; + 1) - 4a; + 2 4x 2 .g(a;) 
/ ( x ) = T^—i — , ,x , -,19 - 24a; ~~ 



[2a;lna; - 2xiP{x + 1) + l] 2 [2a;lna; - 2x?/>(a; + 1) + l] 2 ' 

where; 



I I ^ 2 

g(x) = i)\x + 1) - _ + _ - 24a: 
a; 2a; z 

In [2, Theorem 8], the functions 



tp(x + 1) — In x — 

2a? 



FnW = InrW - - -J In, + * - - ln(2 7 r) - g ^JZ^Fi 

and 

/ x n j 2n+l 

G„(x) = -lnT(x) + f x- - J lna;-a; + -ln(27r) + ^ 



(23) 



(24) 



(25) 



l 2j(2j - l)a; 2 ^i 

for n > were proved to be completely monotonic on (0,oo). This generalizes [16, 
Theorem 1] which states that the functions F n (x) and G n {x) are convex on (0, oo). 
The complete monotonicity of F n (x) and G n (x) was proved in [12, Theorem 2] once 
again. In particular, the functions 



F 2 {x) = In T(x) - [x - i j In x + x - i ln(27r) 

111 1 

+ 



(26) 



and 



12a; 360x 3 1260a; 5 1680a; 7 



Gi(x) = - In r(a;) + [x - ^ J In x - x + ^ ln(27r) 
1 1 1 



(27) 



12a; 360a; 3 1260a; 5 
are completely monotonic on (0,oo). Therefore, we have 

1260a; 5 + 210a; 4 - 21a; 2 + 10 
lnX 2520^ < ^ 



2520a; 7 + 420a; 6 - 42a; 4 + 20a; 2 - 21 

<lna; — — 5 (28) 

5040a; 8 v ' 

and 

210a; 8 + 105a; 7 + 35a; 6 - 7a; 4 + 5a; 2 - 7 , ., . 

210^ < * <*> 

210x 6 + 105a; 5 + 35a; 4 - 7a; 2 + 5 
< 210^ (29) 

on (0,oo). From this, it follows that 
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In x + ip(x + 1) = In a; — ^P{x) 

2x 2x 

1260a; 5 + 210a; 4 - 21x 2 + 10 1 _ 10 - 21a; 2 + 210a; 4 
< 2520a; 6 2x ~ 2520a; 6 ^ ' 

and 

,2 



.,,,.,11 1 (10 - 21a; 2 + 210a; 4 ) 



> 



264600a; 11 

210a; 8 + 105a; 7 + 35a; 6 - 7x 4 + 5a; 2 - 7 1 1 
210a; 9 a? ~ x 

1 (10- 21a; 2 + 210a; 4 ) 2 
+ 2x^ 264600a; 11 
_ 1659a; 4 - 8400a; 2 - 100 
~~ 264600a; 11 

1659(a; - 3) 4 + 19908(a; - 3) 3 + 81186(x - 3) 2 + 128772(a; - 3) + 58679 
~ 264600a; 11 ' 

Hence, the function g{x) is positive on [3, oo). So the derivative f'(x) > on [3, oo), 
that is, the function f(x) is strictly increasing on [3, oo). 
It is easy to obtain 

/fl)-2£^Z> =0.9507..., 

4(487+^1.2-61) = 
•' v ' 5-4 7 -41n2 

3(108 7+ 1081n3-181) 
M ' 5-3 7 -31n3 

This means that the sequence f(n) for n € N is strictly increasing. 
Employing the inequality (30) yields 

, 2520a; 6 ,„ , 12a; 2 (21x 2 - 10) 6 

f(x) > 12a; = - — — > — 

J w 10 - 21a; 2 + 210a; 4 10 - 21a; 2 + 210a; 4 5 

as x — » oo. Utilizing the right-hand side inequality in (28) leads to 
I{X) = ]nx-l/2x-i>(x) ~ 12X2 

1 9 

<- 12a; 

(2520a; 7 + 420a; 6 - 42a; 4 + 20a; 2 - 21)/5040a; 8 - l/2x 

12a; 2 (42a; 4 -20a; 2 + 2l) 

~ 420a; 6 - 42a; 4 + 20a; 2 - 21 
6 

^5 

as x — > oo. As a result, it follows that lim^-nx, f(x) = |. Therefore, it is derived 
that /(l) < f(n) < | for n e N, equivalently, 

2 ( 7 - 12 7) < 1 12n2< 6 

2 7 -l ~ lnn+ l/2n- ip{n + 1) 5 
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which can be rearranged as 

12n» + 2(7 - 12 7 )/(2 7 - 1) ^ ln " + Yn ^ + 1} > 12n' + 6/5 ' 
Combining this with (2) yields (21). The proof of Theorem 1 is proved. 
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